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ABSTRACT 

We consider the one-dimensional random field Ising model, where the spin-spin coupling, J , is 
ferromagnetic and the external field is chosen to be -\-h with probability p and — h with probability 
1 — p. At zero temperature, we calculate an exact expression for the correlation length of the quenched 
average of the correlation function {sQSn) ~ (sq) (Sn) in the case that 2J /h is not an integer. The 
result is a discontinuous function of 2J /h. When p = 5, we also place a bound on the correlation 
length of the quenched average of the correlation function (sQSn)- 
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I. INTRODUCTION 

The one-dimensional random field Ising model is an intriguing example of a system with non-zero 
entropy at zero temperature. The Hamiltonian for this system is 

H ^ -jJ2siSi+i-J2hiSi (1.1) 

i i 

where the spin at each site, Sj, takes the values ihl and > 0. The external field {hi} is frozen 
with the value of the field at each site chosen as an independent random variable with, say, fifty percent 
probabihty to be +/l and with fifty percent probability to be —h. For a fixed external field, {hi], 
and h <: 2J, there is more than one spin configuration, {si}, which minimizes the energy. The zero 
temperature entropy has been calculated^'^'^'^^, and for 2J/h not equal to an integer it depends only 
on the integer q, defined by 

2J , , 

?<-^<? + l ■ (1-2) 

For 2J/h equal to an integer, the entropy is larger than it is when 2J/h is slightly less than or greater 
than this integer. Thus, the entropy is a discontinuous function of 2J//i. For h > 2J, i.e. q — 0, the 
entropy is zero since the spins must follow the external field. 

We are interested in the correlations between the spin at site j and the spin at site j + n. For a fixed 
external field configuration, even at zero temperature, we must average over the different degenerate spin 
configurations. We denote this thermal average by ( ). An object such as {sjSj-^ji) will depend on 
the particular external field configuration, most sensitively on the values of the external field near and 
between the sites j and j -\-n. What is typically measured in scattering experiments is the average over 
the sample 

G{n) = lim — {sjSj+n) (1-3) 

iV— »-oo iV 

3 

where the sum is over the N sites of the system. By the usual ergodic arguments we can replace the 
spatial average by an average over the various possible external field configurations. This allows us to 
write 

G{n)=J^ (1.4) 

where the overbar means average over all external field configurations generated with the probabilistic 
rule introduced earlier for {hi}. 

The ordinary (non-random) one- dimensional Ising model has spontaneous magnetization at zero 
temperature. Any non-zero external random field, chosen with equal probabilities to be ±/l, destroys 
this magnetization. This means that for h ^ 0, 

lim E (sj) = (1.5) 

3 



or alternatively that (sj) = 0. However, for a fixed external field configuration (sj) will generally not 
be zero at the site j. In fact, {sj) and will be correlated so the calculation of 

X(n) = lim ^ [(sjSj+n) - {sj) {sj^n)] (1-6) 

iV— >oo iV 

3 

which is the same as 
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X{n) = (sQSn) - (so) (Sn) (1.7) 

is rather different than the calculation of ^(n)^. 

In this paper we calculate the correlation length, L, of x(7l) which depends on its large n behavior, 
i.e. x(7l) ~ e""'/^. If the external random field is chosen with probability p to he h and with 
probability 1 — p to be —/l we find that 



In 



where q is defined by (1.2) when 2J/h is not an integer. 

We study the correlation function G{n) in the case when the external field is chosen with equal 
probability to be ±/i. We are able to show that for any y such that y > with L given by 1.8 
for p = ^) 'w^s have |G(n)| < for large enough n. We argue, but do not prove, that the bound is 
saturated which would imply that the two correlation functions have the same correlation lengths. 

For both s-nd G(n) we also discuss how our calculations are modified if 2J/h is an integer. 

The random field Ising model is closely related to the random bond Ising model in a uniform field. 
To see this, replace each Si in (1.1) with hiSij h. The new Hamiltonian has a uniform field and random 
bonds Ji,i+i = Jhih%-\-i/h'^ . If the hi are selected with equal probabilities to be ±/i, then the bonds 
take the values ± J with equal probabilities and additionally each bond is independent of the others. 
However if hi — +h with p ^ ^, then the associated random bond model does not have the bonds 
chosen independently. The random bond model with independently chosen bonds is what is usually 
considered in the literature. For example, in references [1,2,3] the entropy is actually calculated in a 
random bond model so these results only apply to the random field model if p = ^ . 

The description of the degenerate configurations which contribute to the zero temperature entropy, 
which we give in section II, is needed before we can attempt to calculate the correlation functions (sections 
III and IV). 

II. SPIN CONFIGURATIONS AT ZERO TEMPERATURE 

For a given external field configuration, {/ij}, the spin configurations, {Sj}, are those that minimize 
the energy (1.1). The first term in the energy favors agreement between adjacent sites whereas the second 
term prefers the spin at a site to agree with the random field at that site. A given {/^i} will not uniquely 
specify the {sj}. However, there can be stretches of the {hi} that force the associated spins to take 
unique values. For example, if we find a very long stretch where all of the hi — +1, then clearly in that 
stretch all of the Si — +1. 

Consider a sequence of k sites at which hi — —1 and imagine that this sequence is flanked on both 
sides by very long stretches where hi = +1. In the flanking regions all Sj = +1. A simple calculation 
shows that if k > 2J/h, then the spins in the sequence all match the external field i.e. Si = — 1, 
whereas if A; < 2J/h the spins are all = +1. Note that if 2J/h is an integer and k is equal to 
2J j h, then these two spin configurations both minimize the energy. Here we see a source of entropy 
present only when 2J /h is an integer. 

To understand how entropy arises for any non-integer value of 2J /h > 1, consider a long stretch 
where hi — +1, followed by a long stretch where hi — —1. This is a deterministic situation where the 
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spins follow the external field. Now at the break point imagine inserting two additional sites where the 
external field takes the values —1 and +1 so the {hi} configuration is 



{hi} = •••+ + + + + - + 

The two spin configurations 

{si} = • • •+ + + + + + + 

and 

{si} = •••+ + + + + 



have the same energy (by symmetry) and by comparing with two other configurations (where the spins 
arc — h and H — in the middle) we see that the two illustrated configurations minimize the energy. This 
is an example of zero-temperature entropy. 

We have seen that regions of constant hi longer than 2J/h force the spins to line up with the field 
and that there are configurations which do not determine the spins. We now state the general rules 
which dictate which regions of the {hi} configuration necessarily determine the spins. 

We denote a region of sites as [£, r] if the left-most site is £ and the right-most site is r. We define 



r 



W[£,r] = ^Y.hi (2.1) 

which measures the difference between the number of sites at which the random field is positive and the 
number at which it is negative in the region [i, r]. We further call [£, r] an region if it meets the 
following three conditions: 



conditions: < 



f(i) W[£,r]> W[£,i] £<i<r 

(ii) W[£,r]> W[i,r] £<i<r (2.2) 
,(iii) W[£,r]>2J/h 



Condition (i) says that starting from £, the number of sites at which hi — +1 minus the number at 
which hi = —1 has a maximum in [£, r] at r and condition (iii) tells us that this maximum exceeds 
(or equals) 2J/h. Similarly, we call [£, r] an R~ region if 



R conditions: < 



'(i) W[£,r] <W[£,i] £<i<r 

(ii) W[£, r] < W[i, r] £<i<r (2.3) 
, (iii) W[£,r] < -2J/h . 



Now an R~^ region favors having all of the spins in [£, r] be +1 over having them all be —1, whereas 
an R~ region prefers all —1 spins over all +1. To guarantee that all spins in an R~^ region [£, r] be +1, 
there should be no R~ subrcgions of [£, r]. We call the region [£, r] a region if it is an R^ region 
with no R subregion. Similarly, we call the region [£, r] a D region if it is an R region with no 
R'^ subregion. 

The spin at a given site will be +1 if the site is in a D+ region. Similarly, if a site is in a D region 
the spin at that site will be —1. Now a given site will be either in a D'^ region, a D~ region, or in a 
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region where the spin is not forced, an E region. With our definitions a D'^ region can be contained 
in a larger D region. We call a D+ region (for maximal) if it is not contained in any other D+ 
region. Similarly, a region is a. D~ region which is not a subset of a larger D~ region. Every 
lattice site is either in a D^, or E region. We cannot have consecutive regions since together 
they would form a D+ region which contained them both. Similarly, two consecutive E regions will be 
considered as one E region. 

To understand what an E region looks like, consider three consecutive regions D+, E and D- 
and let £ be the left-most site of E and r be the right-most site i.e. E — [£, r]. Consider — 1, i] 
as a function of i. Then W[£ — 1,£ — 1] = +1 since £ — 1 is the right-most site of which must 
end in a +1 site. For any k G E, W[£ — 1, A"] < 1 because if this were not the case then the 
region could be extended. Similarly W[r -f 1, r -|- 1] = — 1 and W[k, r -|- 1] > — 1 for /c G E. These 
two inequalities imply that Vl^[£, r] = which means that the entropy region has the same number of 
hi — +1 sites as hi — —1 sites. The function W[£ — 1, i] is equal to 1 at i = £ — 1 and ati — r. It 
can achieve the value 1 but not exceed it at other sites in E, and it also can never go below 1 - 2J/h. 
If it did subregions of the E region would meet the conditions for being D'^ or D~ regions. These 
properties of \V[£ — l. i] will be used when we calculate the correlation functions. 

We now describe the degenerate spin configurations associated with a D^^ED^ region. The spins 
are all +1 in and continue to be -|-1 until some point in E where they switch to — 1 and remain 
— 1 through D^. The last site at which Sj takes the value +1 must be at i = £ — 1 or r or any other 
possible site in E at which W[£ — happens to be +1. We will illustrate this with an example 
momentarily. First note that a mirror construction is used for a D^ED^ region. It is also possible to 
show that no D^ED^ or D^ED^ regions can exist (when 2J/h is not an integer). 

As an illustration, suppose 2 < 2J/h < 3 and we have the {hi} configuration 

{hi} = ••• + + + + + + -+ + ••• 

• •• 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 •• • 

where the numbers below are the site labels. It is useful to plot VF[1, i] from which we can infer the 
values of W on subregions. 

We can see that [1,4] is a D'^ region whereas [11, 15] which has VF[11, 15] = —3 is a D~ 
region. The region [5, 10] is an entropy region and if we look at H^[4, i] for 4 < i < 10 we see that it 
is equal to 1 at i = 4, 8 and 10. The three degenerate spin configurations are 

{s,} = ...+-|--|--|- 

{s,} = ■■■ + + + + + + + + ••• 

{si} = ••■+ + + + + + + + ++ ■■• 

• •• 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 •• • 

It is also interesting to study the same {hi} configuration if 1 < 2J/h < 2. In this case the only 
E region is [9, 10] and the two possible spin configurations are 

{si} = ... + + + + -- + + 

{si} = •••+ + + + -- + + ++ ••• 

• •• 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 •• • 
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which do not coincide with any of the three possibihties for 2 < 2J /h < 3. If 2J /h = 2 then all 
five configurations are degenerate and the entropy is larger than it is on either side of 2J/h — 2. In 
general when 2J/h is an integer there are even more degenerate configurations than those one would 
discover by looking at 2J/h just above and just below its integer value. This is because there can be 
degenerate configurations which within a single E region look in part like those for 2J/h just above its 
integer value and in part look like those for 2J/h just below. 

III. THE CORRELATION LENGTH OF xN 

Recall that 

X{n) = (sosn) - (so) {Sn) (3.1) 

where ( ) is the average over different degenerate spin configurations for fixed {hi} and is the 
average over {hi}. In this section we determine the dominant large n behavior of x(^) for Q < 2J/h < 
q +1 with q an integer. Note that for a given {hi}, Sq and Sfi are either determined by the external 
field or they are not. We can think of Sq and Sn as random variables and {sQSn) — (sq) (Sn) as 
their statistical covariance. If either Sq or is forced by the {hi} to take a particular value then 
i^OSn) ~ {sq) (Sn) vanishes. Thus for {sQSn) — (Sn) (sq) to be non-vanishing, both Sq and Sn must 
be in E regions. However, (sqS^) — (sq) (Sn) also vanishes if Sq and S^ are independent. Now if Sq 
and Sfi are in diff'erent E regions, that is, E regions separated by at least one or Z),^^ region, then 
the value of Sq is independent of the value of Sn and the correlation vanishes. For (sQSn) — (sq) (Sn) 
to be non-zero for a given {hi}, both Sq and Sn must be in the same E region. 

For a given value of n we will calculate the probability, i.e. the fraction of configurations {hi}, 
which have and n in the same E region. The distribution of {hi} configurations is given by assuming 
for simplicity that at site i, hi = +1 ov hi — —1 each with probability one- half. (The calculation is 
easily carried through with ^ replaced by p.) As a first step we calculate the fraction of E regions which 
have length R. Consider an E region which, for example, begins at site 1 and ends at site R and has a 

region to the left and a region to the right. (Note that the sites 1 and R have nothing to do 
with the sites and n mentioned above.) The function Vl^[0, i] as discussed in the previous section has 
the following properties 

(i) iy[0,0] = l ; iy[0,l]=0 ; W[0,R]^1 ; W[0,i]<l for ie[l,R] 

and since W[0, i] > 1 — 2 J/ h with q < 2 J/ h < q + 1 we also have 

(ii) W[0,i] > -q for i e [0,R] . 

The sites just to the right of R form the beginning of a region. Therefore, the function W[0, i] 
fov i > R must take the value —q before it takes the value +1 or else the E region could have been 
extended beyond R. So 

(iii) W^[0, i] iov i > R goes through —q before it goes through +1 . 

The randomly generated field {hi} can be thought of as determining (or as equivalent to) a random 
walk, RW, where position at time i changes by hi/h. The function VF[0, i] with VF[0, 0] = 1 is the 
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position of the random walk at time i given that a,ti — the walk is at + 1 . If we call fji the normalized 
probability that an E region has length R we see from (i), (ii) and (iii) above that 



fjl — N X Prob (RW goes from 1 to 1 in i? steps without hitting +2 or —q) 

(3.2) 

X Prob (RW starting at 1 goes to —q before returning to 1) 

where is the normalization factor. Actually, we will calculate the transform 

oo 

7(A) = E /i^A^ (3-3) 

R=0 

which is more useful for our purposes and from which we can infer /jj. (The R = term in the sum 
corresponds to an E region of zero size which occurs when there are no sites between a region and 
a region.) Note that the transform (3.3) is intimately connected to the transform of the correlation 
function x(A) = S^^=o x{^)^^'j however we only need to calculate (3.3) to infer the large n behavior 
of x(n). 

Turning to the first term in (3.2), let 

Zi(i?) = Prob (RW goes from 1 to 1 in i? steps without hitting +2 or —q) . (3-4) 

To find this, we solve for the more general function 

Zj{R) = Prob (RW goes from j to 1 in i? steps without hitting +2 or —q) (3-5) 

and then set j = 1. Since a walk starts at j and immediately goes to j + 1 or j — 1 we have for 
-q + l<J<l, 

Zm = \Zj-i{R - 1) + \Zj+i{R - 1) . (3.6) 

By (3.5) Z_q{R) = Z2{R) = for i? > and Zi(0) = 1. If we define Z2{-l) = 2 and 
Zj{-1) = for -q<j< 1, then (3.6) holds for = as well as for R > 0. 
We define the transform 

oo 

Zj{X)= E Zj{R)X^ (3.7) 

R=-l 

which from (3.6) gives for —q + 1 < J < 1, 

Zj(X) = ^AZ,_i(A) + lxZj+,(X) . (3.8) 
We can solve (3.8) by making the ansatz that 

Zj(X)^au^ + ^v^ (3.9) 
with the boundary condition that Z_^(A) — and .^2 (A) = 2/ A as explained above. The solution is 

2 



Zj{X) ^ 



u 



9+2 _ .,9+2 



(3.10) 
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with 



Thus we obtain 



1 



and 



A 



A2 



^9+1 _ yq+l 



(3.11) 



(3.12) 



_ yq+2 

We now return to (3.2) and we see that the second probabiUty factor is independent of R and can 
therefore be absorbed in the normahzation factor N. Thus we have for the transform /(A) defined by 
(3.3), 

/(A) = NZi (A) (3.13) 



which by (3.12) gives 



/(A) = 



27V /m^+i -yg+i^ 



(3.14) 



A \ui+'^ - vi+'^ J 

Now (3.14) can be expanded in only non- negative powers of A as in (3.3). The normahzation condition 

oo 

J2 Jr = 1 is equivalent to /(I) = 1 which allows us to solve for N and we obtain 
R=0 



/(A) 



q + 2 1 ( u 



g + 1 A \ui+^ -vi+^ 



(3.15) 



Again, if we expand /(A) as a power series in A, the coefficient fji is the probability that an E region 
has length R. 

We now turn to finding 



Qn — Prob(0 and n are in the same E region) . 

■R-n- 



R 



(3.16) 
(3.17) 



The answer is 

= TV' ^ RfR 

R>n 

where we explain each factor in turn. The factor N' contains the probability that is in an E region 
and other n-independent factors. For the E region to contain and n it must have length R > n. The 
factor R fji is proportional to the probability that an E region has length R given that is in it. The 
factor [R — n) /R is the probability that n is in an E region of length R given that is in it. 
We can also define the transform 



which by (3.17) is 



Q(A) = E QnA^ 
n=0 



_ oo R-l 

R=l n=0 



It is straightforward to do the sum on n and then on R to obtain 



A) 



dX 



A 

A=l~ (1^ 



/(I) 



A 

(1^ 



/(A) 



(3.18) 



(3.19) 



(3.20) 
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We have the expUcit form of /(A) through (3.15) so the coefficients in (3.18) can be determined for 

all n. Thus we have computed the (unnormalizc^d) probability that and 71 are in the same E region. 

The large Tl behavior of Q^i can be extracted if wc know the smallest value of A > 0, say A*, at 
which (3.20) blows up. This is because the expansion (3.18) will blow up first at A = A,,, if Qn ~ A"*^. 
Now (3.20) does not blow up at A = 1 as can be seen by expanding /(A) about A = 1. The only way 
for (3.20) to blow up at A 7^ 1 is for /(A) to blow up. Prom (3.15) we see that this can occur only if 
u'^'^^ — v^"*"^ where again u and V are given by (3.11). A simple calculation gives 

from which we conclude that, for large n, the probability that sites and n are in an E region goes as 

Qn ~ cos" (-^] ■ (3-22) 



Given that and n are in the same E region we need to calculate {sQSn) — {sq) {sn) where the 

E on the overbar denotes average only over those {hi} for which and n are in the same E region. By 

examining the degenerate spin configurations in an E region one can see that (so^n) ~ (*o) (*n) ^ 

for all {/ii}, so no cancellations take place; in this average. We are interested in n large and the most 

probable configurations contributing to this average are those for which the E region is just slightly 

longer than n so the site and the site n arc near the edges of E. The number of degenerate spin 

configurations associated with an E region is proportional to its length from which we can estimate that 

E 1 

{sQSn) - {so) (Sn) ~ ^• 

By combining (3.21) with the estimate of the previous paragraph gives, for large n, 

vr 



X{n) = {sQSn) - {so) {Sn) ~ COs" I ^ j ^^'^^^ 



from which we infer that the correlation length is 

L = ^ ^ . (3.23) 

In cos (^j 

It is worth noting that if we add a constant external field, no matter how small, then we destroy 
the zero temperature entropy since the degeneracy is lifted. In this case x(7l) = and there is also a 
non-zero magnetization, i.e. {sq) 0. Alternatively we can pick the random field at each site to be 
-\-h with probability p and to be —h with probability 1 — p, and then (sq) 7^ unless p = ^. In this 
case it is straightforward to redo the calculation of the correlation length and we get (1.8). 

If 2J/h is an integer, say k, then the calculation of the large n behavior of x(7l) changes in two 
ways. First, the probability that sites and n are in the same E region is larger when 2J/h is equal 
to k than when 2J/h is slightly greater than k. For example when 2J/h — 1, site and site n are in 
the same E region in the following configuration: 

{hi} = ... + + - + - + - + - + + ... 

n 
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whereas the entire pictured region is if 1 < 2J/h < 2. But this only approximately doubles the 
chance that and n arc; in the same E region and has no effect on the correlation length. 

However the typical value of (sQ'^ra) ~ (-^o) i^n) a-lso changes if 2J/h = k as opposed to k < 
2J/h < k + 1. In an E region of length n when 2J/h is not an integer there are of order n 
configurations and (so^n) ~ (^o) (Sn) is of order Tl ^. When 2J/h is an integer there are more 
configurations. For example when 2J j h = \ aji E region of length L has Eij^2 configurations where 
Ft, is the L-th Fibonacci numberl^l (Fi = F2 = 1 ; Fij.2 = Fl+1 + Fl)- Now for each {hi}, 

^ 

{sqSh) — (-^o) i^n) is still non- negative but {sqSji) — (-^o) (^n) is of order F~^. For 1 < 2J/h < 2 
we have that x(^) ~ (^) can be seen from (3.22). For 2J/h = 1 we have x(n) ~ (| j F~'^ 

and since Fn ^ ((■'-"'" ^) /^) ^^^^^ ^^^^ L"-*^ = In ((1 + V5) /2J. 
IV. A BOUND ON THE CORRELATION FUNCTION G(n) 

We are interested in the correlation length of {sQSf^) in the case when the external field is chosen 
with equal probabilities to be ±/i at each site. In this case (sq) — 0. We begin by using the symmetry 
of the problem to identify a class of the {/ii} which has the property that (so^n) averaged over this 
class is zero. Roughly, this is the class of {hi} where sites and n are separated by at least 2 disjoint 
D regions. We then will estimate the probability that the sites and n are not in this class. This turns 
out to have the same large n behavior as the probablity that sites and n are in the same E region, 
which was relevant in calculating 

Suppose we are given a particular external field configuration {hi}. Let [Eq, Tq] be a minimal D+ 
or D~ region with Tq ^ and rg as small as possible. (A minimal or D~ region has no subregion 
which is a D'^ or D~ region.) Let [i^i, Tn] be a minimal D'^ or D~ region with ^ ^ 
as large as possible. Suppose [^Oj'^o] and [^td '"n] are disjoint, that is Tq < (which is likely if n 
is large). Switching all the random field signs at in, i-n+l-, in+2, ■ ■ ■ has the effect of leaving (sq) 
unchanged but reversing the sign of (Sn)- Note that (sq) (Sn) — (sQSn) in these cases and also that 
performing the switch twice returns us to the original configuration. Thus (sQSn) averaged over all of 
these configurations is zero. 

Next we estimate the probability of obtaining a configuration of external fields, {hi}, with Tq ^ in- 
It is these configurations which produce a non-zero correlation function. The reader who wishes to skip 
the details of this estimate should proceed to equation 4.11 which gives the result. 

We begin by calculating the distribution of Tq, which is the smallest non- negative site which is 
the right-most site of a D region. Fix h—i, h—2, h—s ■ ■ ■ ■ We actually calculate the distribution 
of To conditional on these values. We will see, however, that the probability of rg being large will be 
essentially the same for any choice /i— 1, /i— 2; ^—3 • • • • Given the random field at the negative sites, 
find the largest value of 0, < so that [b, a] is a Z) region for some h <. a. Without loss of generality 
assume it is a region. (Note that the external field at the non-negative sites may make [6, a] part of 
an even larger D region. However a is defined only using the values of the fields at the negative sites.) 

Consider l^[a + 1, r] as a function of r > a + 1. If Vr[a + 1, r] reaches the value 1 before it goes 
through —q—\, then the region \b, a] can be extended. This cannot happen for r < for if it did 
a would not be the largest negative site ending a D region on the right. If VF[a + 1, r] reaches the value 
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—q — 1 at some r before reaching the value 1 then there is a D region with r as its right end. Again, by 

assumption, this cannot happen for r < 0. We can sec now that Tq is the smallest value of r > such 
that iy[a + l,r] = 1 orW[a + l,r] = -q-1. Because W [a + l,r] = W[a + 1, -1] + W[0,r] 
for r > 0, we can say that Tq is the first r such that W[0, r] = ^4 or W[0, r] = B with 

A^l-W[a + 1,-1] (4.1) 
B ^ -q-1- W[a + 1,-1] 

(where for a = —1 wc define W^[0, —1] = 0). Thus for a fixed configuration at the negative sites, 
the distribution of Tq > depends only on the single number W[a + 1, —1] which obeys —q < 
W[a + 1,-1] < 0. 

As in the previous section we view H^[0, r] as equivalent to a random walk. Let 

Yj{R) = Prob(RW starting at j first hits A > or S < at step R) (4.2) 

Note that 

Prob (ro = I h-i, h-2, ■■■) = Yo{R + 1) . (4.3) 

Now 

Yj{R) = ^Yj-iiR - 1) + ^Yj+iiR - 1) (4.4) 

with the bomidary conditions that 1^(0) = Yb{0) = 1 and Ya(-R) = Yb{R) = for > 0. We 
can solve for the transform of Yj(^R), 

oo 

YjW = E U^)^"" (4-5) 
R=0 

as we did in the previous section to obtain, 

^0^^) = ui+^ - vi+^ ^^-^^ 

where again u and V are given by (3.11) and we have used the fact that A — B — q + 2. Note that 
Yq[X) blows up for the first time at A* given by (3.21) which is independent of A and B. 

The coefficient of A'^+^ in (4.6) gives the probability that Tq has the value of R given a fixed 
h^l, /i— 2 • • • which determine A and B. Similarly wc could obtain an identical expression for the 
probability that Tl — £n has a given value for a fixed hji^\, • • • • Wc arc interested in calciilating 

the probability that Tq > £n with h—i, h—2 . . . = {/i<} and hn+1, hn+2 ■ ■ ■ = {^>} both fixed. 
Now 

Prob(ro > £n I {h<}, {h>}) = 1 - Prob(ro < £n I {h<}, {h>}) 

= 1-5^ Prob (ro ^i,n-in^j\ {/i<}, {/i>}) . (4.7) 

i,j>0 

i+j<n 
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The probability that rQ — i depends on the random field at sites < i while the probability that 
n — in = j depends on the sites ^ n — j which do not overlap in the sum in (4.7) so the distribution 
can be taken as independent. It then follows that 

Prob (ro >in\{h<}, {h>}) = P'^ob (ro = i | {/i<}) Prob {n-in^j\ {h>}) (4.8) 

hj>0 

i+j>n 

so for the purposes of our calculation we can treat the full distributions of ro and n — as independent. 
Consider the transform of the probability that Tq + n — £n has the value k: 

oo / oo \ 

^ Prob {ro + n-£n^k\ {h>}, {/l<}) A'^ = ^ Prob {tq ^ i \ {h<}) A' X 

k=0 \i=0 / 

00 \ 

Y,Proh{n-in^j\{hy})Xn .(4.9) 
\j=o J 

Both transforms on the right hand side are of the form (4.6) and the product blows up at A* given by 
(3.21) so we can say that for k large 

Prob(ro + n-4 = A; I {h<},{h>}) ^ cos'' ( (4.10) 



from which we infer that 



for large n. 



Prob (ro > in) ~ cos" f j (4-11) 



The only configurations of the external field which contribute to {sQSn) arc those for which Tq ^ 
Given a configuration with rg ^ in we expect (sQSn) ~ +1 or — 1. We have not shown that 
cancellations do not conspire to make the average of {sQSn) over those configuration with ro > in of 
order with \x\ < 1. Hence we can only assert that 



TT 



if y > cos , then {sQSn) < , for n large enough. (4-12) 



We can write (sQSn) as 



{sQSn) = x(n) + (so) {Sn) , (4.13) 

we know that x(^) ^ 0> and from (3.23) we see that it decays as COs" (^^) • Therefore if (sq) {sn) 
is non-negative we can conclude that (sQSn) has as its correlation length, L, given by (3.23). However 
we have not been able to prove that {sQSn) ^ although the following argument makes us believe 
that it is. Consider setting the spin-spin coupling, J, equal to zero, which gives {sQSn) = 0. For J > 
we expect the feromagnctic coupling to induce a positive correlation between Sq and Sn, even in the 
quenched average. For this reason we believe, but have not proven, that the correlation length of {sQSn) 
is L given by (3.23). 

Finally we remark that if 2J/h — k, an integer, the argument leading to (4.11) again gives (4.12) 
with q — k. 



-11- 



ACKNOWLEDGEMENTS 

We thank R. J. Birgeneau for tempting us with this problem and A. N. Berker and M. Kardar for 
useful discussions. 

REFERENCES 

1. M. Puma and J. F. Fernandez, Phys. Rev. B 18 (1978) 1391. 

2. A. Vilenkin, Phys. Rev. B 18 (1978) 1474. 

3. B. Derrida, J. Vannimenus and Y. Pomeau, J. Phys. C: Solid State Phys. 11 (1978) 4749. 

4. B. Derrida and H. Hilhorst, J. Phys. C: SoUd State Phys. 14 (1981) L539. 

5. J. K. Williams, J. Phys. C: Solid State Phys. 14 (1981) 4095. 



-12- 



